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The S = 1/2 Heisenberg chains with bond alternation and randomness on the strong bonds 
are studied by the density matrix renormalization group method. It is assumed that the odd- 
th bond is antiferromagnetic with strength J and even-th bond can take the values Ja and 
Jf (Ja > J > > Jp) randomly. The ground state of this model interpolates between the 
Haldane and dimer phases via a randomness dominated intermediate phase. Based on the 
scaling of the low energy spectrum and mean field treatment of the interchain coupling, it 
is found that the magnetic long range order is induced by randomness in the intermediate 
regime. In the magnetization curves, there appears a plateau at the fractional value of the 
saturated magnetization. The fine structures of the magnetization curves and low energy 
spectrum are understood based on the cluster picture. The relation with the recent experiment 
for (CH3)2CHNH3Cu(Cl:^Bri„^)3 is discussed. 
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1. Introduction 

In the recent studies of quantum many body problem, 
, the one-dimensional random quantum spin systems have 
' been attracting a renewed interest from theoretical and 
, experimental viewpoints. ^^^^^ 

Among them, the phenomenon of the disorder induced 
, order have been widely investigated. Experimentally, 
' Uchinokura and coworkers'^' have found the antiferro- 
, magnetic ordered phase in Zn, Mg, Si-doped spin-Peierls 

■ compound CuGe04 as one of the earliest examples of 
, this type of phenomenon. The theoretical explanation is 
' given by Fukuyama and coworkers^^ using the bosoniza- 

tion approach. 

Similar phenomena are observed in Zn-doped SrCu203^^ 
1 which is the quasi- 1-dimensional S = 1/2 ladder sys- 

■ tem and Mg-doped PbNi2V208*'^'' which is the quasi- 
, 1-dimensional = 1 Haldane gap system. Correspond- 
' ingly, the effect of the bond and site randomness on the 
, spin gapped quasi- 1-dimensional S = 1/2 and S = 1 an- 

■ tiferromagnets are studied by the quantum Monte Carlo 
, method by Yasuda and coworkers.*^ They have also 
' found the randomness induced long range order for ap- 
, propriate range of randomness. 

Recently, Manaka and coworkers^' studied the 
magnetic and thermal properties of the quasi-one- 
dimensional compound (CH3)2CHNH3Cu(Cl2;Bri_2;)3 
(IPACu(Cla;Bri_a;)). For a; = 0, this material is 
the S = 1/2 antiferromagnetic-antiferromagnetic alter- 
nating Heisenberg chain (AF-AF chain) whose ground 
state is the dimer phase^^' while for a; = 1 it is 
the 5 = 1/2 ferromagnetic-dominant ferromagnetic- 
antiferromagnetic alternating Heisenberg chain (F-AF 
chain) with Haldane ground state. ^^' In the interme- 
diate concentration regime, however, both the Haldane 
and dimer phases are destroyed by randomness. Re- 
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markably, Manaka and coworkers^-* experimentally found 
the magnetically ordered phase for 0.44 < a; < 0.87. It 
should be also remarked that the energy gap estimated 
from the temperature dependence of the specific heat 
and susceptibilty remains finite even in the close neigh- 
bourhood of the critical concentration. 

Modivated by this experiment, we investigate the ran- 
dom S = 1/2 Heisenberg chain whose Hamiltonian is 
given by, 

N N 

H — ^'^^ JS2i-lS2i +''^^ JiS2iS2i+l, (1) 

1=1 1=1 

where J > and Ji — Jf(< 0) with probability 1 — p 
and Ji = Ja(> 0) with probability p. The ground state 
of this model interpolates between the Haldane phase 
{p = 0) and dimer phase (p — 1). 

As discussed by Manaka and coworkers, the rela- 
tion between x and p is not trivial because the exchange 

paths between two Cu ions are bibridged bonds. For 

/-II 

a; = 1, the ferromagnetic bonds are Cu <q> Cu bonds, 
and for x = 0, the strong antiferromagnetic bonds are 

Br 

Cu <T3 > Cu bonds. The addition of the Br ions into 
Br 

Br 

the a; = 1 chain induces the Cu <qJ> Cu bonds which 
are absent in the pure chains. We assume, however, 
that the Cu <q> Cu bonds are strongly antiferromag- 
netic because otherwise the Haldane phase cannot be 
destroyed and no magnetic order can take place for rela- 
tively small Br concentration (1 — a; ~ 0.13) as observed 
in the experiment. In this case, the concentration of 
the ferromagnetic bond 1 — p is equal to x^. Further, 
we may safely ignore the Cu <gj.> Cu bonds in the 
concentration regime 1 — x << 1 which we are mainly 
concerned in this paper. Actually at the critical con- 

Br 

centration x = 0.87, the concentration of Cu <gj.> Cu 
bonds is (1 — x)"^ ~ 0.017 which is sufficiently small. We 
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can therefore set p = 1 — in this regime. In the fol- 
lowing, we take Jf = — 2J following the experimental 
estimation^) and J = 1 to fix the energy unit. 

In the present system, the spin gap state in the pure 
system is intrinsic. Therefore the mechanism of the dis- 
order induced order is not related with the lattice de- 
grees of freedom as in the spin-Peierls systems.'^"''') The 
situation is somewhat similar to the site depleted 2- 
dimensional 5=1 antiferromagnet studied by Yasuda 
and coworkers^) if Jp and Ja bonds are infinitely strong. 
However, as explained later, the finiteness of these bonds 
introduces interesting fine structures in the single chain 
properties which should be observable in experiments. 
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Fig. 1. The system size-dependence of < In A > for (a) Ja = 2 
and (b) Ja = 4 with Jp = —2. The solid line is the fit to the 
relation In A = const. — zlnN. 

This paper is organized as follows. In the next sec- 
tion, the single chain properties of the present model are 
discussed. In subsection 2.1, the low energy spectrum of 
the single chain is calculated and shown to have quantum 
Griffiths singularity. The physical origin of this behav- 
ior is also explained in the cluster picture. In subsection 
2.2, the magnetization curve of this model is presented 
and the randomness induced plateau is shown to appear. 
The fine structures of the magnetization curve are also 
explained in terms of the cluster picture. In the third sec- 
tion, the effect of interchain coupling is studied by means 



of the interchain mean field approximation. It is shown 
that the disorder induced magnetic order takes place in 
appropriate concentration regime. The final section is 
devoted to summary and discussion. The brief report on 
the main results of this work is published in ref. 11. 

2. Single Chain Properties 

2.1 Low Energy Excitation Spectrum 

To support the long range order in the coupled chain 
systems, enough number of low energy states are required 
in the spectrum of the constituent single chains. We 
therefore calculate the size dependence of the energy gap 
A of the random ensemble of spin chains described by 
Eq. (1). Due to computational reason, we concentrate 
on the regime of small p (small 1 — a;) in the following. 
Figure 1 shows the system size dependence of the average 
of the logarithm of the energy gaps A with Jp = — 2 and 
(a) Ja = 2 and (b) 4 for 10 < 2N < 80 calculated 
by the density matrix renormalization group (DMRG) 
method. -^^'-^^^ The average is taken over 400 samples. 
In this parameter range, these curves are fitted well by 
< In A >~ const. — ^InA'' where z defines the dynamical 
exponent. This implies that the energy gap is scaled by 
which is a typical behavior of the quantum Griffiths 
phase. ^'^"^ If we assume such scaling holds for all low 
energy excitation spectrum, the distribution function of 
the low lying excitation energies scales as 

P{A)dA = N''f{N'A)dA (2) 

where P{A)dA is the number of states with excitation 
energy in the range [A, A -|- dA]. For large N, P{A) 
should be proportional to the system size N, so that 
f{x) as a; ^ oo. Thus, in the thermodynamic 

limit, we have 

P{A)dA oc TVA^/^-^dA. (3) 

Therefore the number of the low energy states diverges 

if z > 1. As shown in Fig. 2, z increases with p and 
becomes larger than unity above a critical value Pc where 
Pc 0.3 for Ja = 2 and Pc — 0.1 for Ja = 4. We can 
expect that the long range order would be stabilized for 
p > Pc if the interchain coupling is switched on. 
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Fig. 2. The p-dependence of the dynamical exponent z. The 
circles axe the estimation from the system size dependence of 
< In A >. The lines are z = aj \ Inp | where a is determined 
from the q dependence of Aq. 

These features of the low energy spectrum can be un- 
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Fig. 3. The g-dependence of Aq. The inset shows the g-cluster. 



derstood in the following way. Let us consider a cluster 
consisting of q JA-bonds and q — I J-bonds embedded 
in the F-AF chain as depicted in the inset of Fig. 3. 
This is called 'g-cluster' in the following. The middle 2q 
spins form a tightly bound singlet cluster. The two spins 
connected to both ends of this cluster by J-bonds are al- 
most free but weakly coupled with each other mediated 
by the quantum fluctuation within the strongly coupled 
cluster. For small p, we can regard the whole system as 
a random assembly of g-clusters. The g-dependence of 
the singlet triplet gaps A, of g-clusters calculated by the 
DMRG method is well fitted by A, ~ Aoe""i as shown 
in Fig. 3. On the other hand, the number of g-clusters 
in a chain is proportional to Np'^{l — p)^. Eliminating 
q, the number P{A)dA of the g-clusters with energy gap 
between A and A + dA behaves as Eq. (3) with 



z = a/ \ Inp I . 



(4) 



This value of z is plotted by the solid and dotted lines 
in Fig. 2. The agreement with the values obtained by 
fitting the numerical data of energy gap in Fig. 1 is 
good for small p. As p becomes larger, the interference 
between the clusters would prevent the above simple- 
minded interpretation. 

Using the formula (4), we can estimate the parameter 
regime in which the low energy spectrum has divergent 
singularity as p > Pc = exp(— a). The critical value Pc 
is plotted against Ja in Fig. 4 using the values of a 
obtained from Fig. 3. For p > pc, it is possible that the 
long range order is stabilized in the presence of interchain 
coupling. 

The above picture explicitly demonstrates that the low 
energy excitation of the present model is dominated by 
the large size clusters and the characteristic size of the 
clusters increases by the power law with the decrease of 
the energy scale. This is a typical feature of the Griffiths 
phase. This singular excitation spectrum is reflected in 
the low temperature magnetic susceptibility x and mag- 
netic specific heat C as x ~ T^/^-i and C ~ T^^. 




Fig. 4. The parameter regime in which the low energy spectrum 
has divergent singularity. 
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Fig. 5. Magnetization curves for (a) Ja = 4 and (b) Ja = 2 with 
Jp = —2 and 2N = 100 for various values of p. The magnetic 
field and the magnetization per site are plotted in units of J/s/Ub 
and g/L^B) respectively. 



4 



Kazuo HiDA 



2.2 Magnetization Curve 

The magnetization curves at T = is calculated by 
the finite size DMRG method for 2A^ = 100 for Ja = 2 
and 4 with Jp = —2 as shown in Fig. 5(a) and (b), re- 
spectively. The average is taken over 100 samples. A 
big plateau appears at M = Ms(\ — p) where Mg is the 
saturated magnetization and small jumps appear below 
and above the main plateau. Again, the physical inter- 
pretation of this structure can be given in terms of the 
f^-clustcrs. On the main plateau the spins connected by 
JA-bonds form singlet clusters and remaining spins are 
polarized along the direction of the magnetic field. These 
plateaux are well quantized so that it must be useful to 
determine "p directly from experimental data. Similar 
fractional plateaux has been also found in random poly- 
merized XXZ chains. 

Recently, Totsuka^^^ analytically discussed the effect 
of randomness on the magnetization plateaux. At first 
sight, our fractional plateau appears to contradict with 
his criterion (Eq. (36) of ref. 13). However, his argument 
concerns the fate of the plateaux which already exist in 
the regular system and does not rule out the plateaux 
induced by the bond randomness in our model as well as 
in the model in ref. 12 as commented by Totsuka himself 
at the end of his paper. 

The small jumps above the plateau are the contribu- 
tion from the spins connected by the Ja bond and those 
below the plateau are the contribution from the pairs of 
almost free spins which support the low energy magnetic 
exciations. The position of these jumps are identified by 
the DMRG calculation for the g-clusters as indicated by 
the up and down directed triangles in Fig. 5(a) and (b) 
for each q. The position of the main jumps above the 
plateau _ffaq and those below the plateau iJbq due to the 
g-cluster are depicted in Fig. 6 against Ja for g = 1 and 
2. Because is sensitive to Ja, we could determine Ja 
from experimental value of iJai- However, considering 
that the saturation field of the x = \ compound (p = 0) 
is already above 40T,^^^ the observation of this jump 
would be rather diSicult within the presently available 
pulse magnetic field ~ SOT unless Ja is relatively small. 
Although less precise, Ja could be also determined from 
experimental data for iJbi which should be observable 
within the presently available magnetic field. In re- 
cent magnetization measurement,^^ however, no promi- 
nant structure is observed in real IPACu(Cla;Bri_x) for 
X > 0.87 in the low field regime. The origin of this dis- 
crepancy is unclear at present. 



3. Effect of Three Dimensionality ■ 
duced Magnetic Order 



Disorder In- 



As demonstrated above, the low energy sectors of the 
present model is dominated by the almost free spins in 
g-clusters with large q. In the presence of the interchain 
coupling, we therefore expect that the 3 dimensional net- 
work of these spins sustain the magnetic order observed 
in the experiment mediated by the interaction with spins 
in neighbouring chains. Let us assume the Hamiltonian 




Fig. 6. The JA-dcpcndcncc of the magnetic fields where the jumps 
due to (jr = 1,2 clusters take place (a) above the plateau (solid 
line: -ffali two dotted lines: -ffa2) B,nd (b) below the plateau 
(solid line: -ffbl) dotted line: -ffb2)- ^^F = —2. 



with the interchain coupling as follows, 

N N 
j i=l 



i=l 



2N 



(5) 



»=i <j,j'> 

where j and j' distinguish the chains and J2<j j'> 
notes the summation over nearest neighbour chains. We 
assume that the correlation between two almost free 
spins separated by the Ja bonds is mainly mediated 
by the F-AF part of the neighbouring chains, because 
the correlation length of the AF-AF part is small com- 
pared to that of the F-AP part. Furthermore, for small 
the probability to find the Ja bonds on the neighbour- 
ing chains between two almost free spins is small. As 
a result, the interchain interaction and the short range 
correlation within the neighbouring chains support the 
ttii-type long range order as far as p is small. 

For quasi-one-dimensional systems, the mean field ap- 
proximation for the interchain coupling gives reliable re- 
sults, because a large number of spins are envolved 
in the low energy long wave length fluctuations in each 
chain. These fluctuations are frozen by interchain inter- 
action resulting in the three dimensional ordering. We 
therefore employ the interchain mean field approxima- 
tion assuming the TT].]. order to obtain the single chain 
mean field Hamiltonian given by. 



N 



N 



2N 



jjMF _ ''^JS2i-iS2i + ''^JiS2iS2i+l — (6) 



i=l 



i=l 



i=l 



with hi = —ZcJ± < SI > and Zc is the number of near- 
est neighbour chains. Reflecting the ttii order, we take 
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h4i = /i4j+i = h and /i4i+2 = /i4i+3 = —h. In Fig. 7, the 
staggered magnetization nist =< jw^i^ii^^yi^ii + 
^ii+i) plotted against A = ZcJ± for (a) Ja = 2 and 
(b) Ja = 4 with Jf = -2. The chain length is 27V = 200 
and finite size DMRG method is used. The average is 
taken over 200 samples. For p — 0, nisi vanishes unless 
A is larger than a critical value Ac. For finite p, however, 
mst takes a finite value even for small A < Ac. The mag- 
nitude of mst increases with p and Ja) namely with the 
increase of randomness. It should be noted that the long 
range order for A < Ac starts to appear around p ~ 0.3 
for Ja = 2 and around p ~ 0.1 for Ja = 4. Theses val- 
ues are approximately consistent with the estimation of 
Pc from the energy gap scaling. Therefore the TTii-type 
long range order is stabilized for appropriate strength of 
the interchain coupling if the low energy density of state 
has divergent singularity as expected. 
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Fig. 7. The A-dependence of staggered magnetization for (a) Ja = 
2 and (b) Ja = 4 with Jp = —2. The solid line is the staggered 
magnetization for p = 0. 

In the experiment, the long range order is observed for 
X < 0.87. If we assume p = 1 — x^, the regime x < 0.87 
corresponds to p>^0.24. Although the exact value of Ja 
is unknown, we find from Fig 7 that the long range order 



is stabilized in this concentration range even for A < Ac 
both for Ja = 2 and 4. These results explain the qual- 
itative features of the experimental observation of the 
magnetic ordered state. For the quantitative compari- 
son with experiments, the strength of the Cu <(^|> Cu 
bonds need to be determined. Similar analysis is made 
for the random 5=1 Heisenberg antiferromagnet by 
Villar et aZ.") 

4. Summary and Discussion 

The S = 1/2 Heisenberg chains with bond alternation 
and randomness on the strong bonds are studied by the 
DMRG method. The low energy spectrum is shown to 
have the Griffiths type singidarity and the physical ori- 
gin of this behavior is explained based on the g-cluster 
picture. In the magnetization curves, there appears a 
randomness induced plateau at the fractional value of 
the saturated magnetization. This plateau and the fine 
sructures of the magnetization curves are also under- 
stood based on the q-cluster picture. By the mean field 
treatment of the interchain coupling, the magnetic long 
range order is shown to be stabilized by randomness 
in the intermediate concentration regime. The results 
are discussed in relation with the recent experiment for 
(CH3)2CHNH3Cu(CUBri_,)3. 

We carried out the interchain mean field calculation 
assuming the tTii type long range order. However, the 
possibility of different type of ordering such as spin glass 
ordering cannot be excluded. From this viewpoint, the 
experimental determination of the magnetic structure by 
neutron scattering experiment is also hoped to elucidate 
the nature of the long range order. 

Recently, Nakamura^^^ has shown that the titi corre- 
lation becomes most critical in the intermediate concen- 
tration regime in the present model from the nonequilib- 
rium relaxation method analysis of the quantum Monte 
Carlo data. We have also checked the possibility of the 
TiTi order within the DMRG and interchain mean field 
approximation. According to our calculation, the stag- 
gered magnetization for |||| order is much lower than 
that for the ttii order. For example, for Ja = — Jf = 2 
and p = 0.4, the staggered magnetization for the 
order is approximately one order of magnitude smaller 
than that for the ttii order. For other values of param- 
eters Ja, Jf and p, the ratio of titi order to ttii or- 
der is even smaller. Especially, the titi order decreases 
rapidly with the decrease of p and becomes less than 
10-3 for p < 0.3 with Ja = - Jf = 2 and for p < 0.2 
with Ja = 4, Jp = — 2 within the regime A < Ag. There- 
fore we may conclude that the dominant order is of the 
ttii type although it is possible that the weak titi order 
also coexist. It should be noted that Nakamura's calcu- 
lation does not explicitly include the effect of interchain 
coupling. In the random systems, the criticality of the 
correlation in a single chain does not always imply the 
corresponding long range order in its quasi- 1- dimensional 
counterpart. 

The absence of the fine structure in the low field mag- 
netization curve of the real IPACu(Cl2;Bri_j;) is not im- 
derstood in the present calculation. In the presence of 
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interchain coupling, however, it is possible that the al- 
most free spins form localized singlet clusters mediated 
by interchain coupling. This can suppress the low field 
structure in the magnetization curve. However, such lo- 
cal correlation effect is not properly described within the 
present interchain mean field approximation. This prob- 
lem is left for future studies. 

In this work, we concentrated on the ground state 
properties due to the limitation of the numerical method 
(DMRG). The finite temperature effects must be impor- 
tant for the direct comparison with experiments. This is 
left for future studies. 
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